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Abstract 



We consider topological supercurrent excitations (SC) in ID mesoscopic rings. 

> " 

ON ' Under certain conditions such excitations are well-defined except for (i) a tun- 

o 

O 

' which may be characterized by the amplitude A, and (ii) a decay of SC as- 



neling between resonating states with clockwise and anti-clockwise currents, 



sisted by phonons of the substrate, both effects being macroscopically small. 
Our approach being based on the hydro dynamical action for the phase field 



and its generalization to the effective Hamiltonian explicitly takes into ac- 

o ' 

' count transitions between the states with different topological numbers and 

turns out to be very effective for the calculation of A and estimation of the 
decay width of SC, as well as for the unified description of all known ID 
superfluid-insulator transitions. 

Most attention is paid to the calculation of the macroscopic scaling of A 
(the main superfluid characteristic of a mesoscopic system) under different 
conditions: a commensurate system, a system with single impurity, and a 
disordered system. The results are in a very good agreement with the exact- 
diagonalization spectra of the boson Hubbard models. 

Apart from really ID electron wires we discuss two other important ex- 
perimental systems: the 2D electron gas in the FQHE state and quasi-lD 
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superconducting rings. We suggest some experimental setups for studying 
SC, e.g., via persistent current measurements, resonant electro-magnetic ab- 
sorption or echo signals, and relaxation of the metastable current states. 



I. INTRODUCTION 



The problem we address here is a general one for the condensed matter physics and arises 
each time when we deal with the sample of nontrivial topology. Namely, if we write down 



with the topology of a ring we may define the topological excitations characterized by the 
winding number 



where the integration contour T is around the ring hole. Clearly, these excitations exist in any 
space dimension and for 1^0 carry current due to nonzero gradient of the phase field. The 
stability of topologically excited states does depend on the dimensionality and properties of 
the bulk system. In the normal state of the 3D and 2D liquid the relaxation time of the SC 
due to disorder or crystal potential is extremely short (of order of the transport scattering 
rate). However below the superconducting or superfluid phase transition the relaxation 
mechanism of I is only through the vortices, which are macroscopic objects in a sense that 
their energy is system size dependent (the pairs of vortices have finite energy, but their only 
effect on / is to provide short-time fluctuations in / leaving the long-time average of the 
topological quantum number unmodified). In the ground state or at a very low temperature 
the probability of vortex creation by quantum or thermal fluctuations can be ignored, and 
the integral (|1]) is a true invariant of the system. Note also, that at finite temperatures 
below the critical point there may be introduced a regularized operator of phase in terms of 
which / is also an exact (with macroscopic accuracy) constant of motion [I]]. 

The above argument does not work in ID, because one can not define vortex in one 
spatial dimension. It means that the dynamics of / is defined by microscopic parameters 
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the field operator in terms of density and phase \l/(x) =| p(x) 



/2 e «$(x)^ thgjj j n a sam ple 
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of the problem, and the bare transition amplitude between different quantum numbers is 
independent of the system size. Calculating the relaxation of the topological excitations 
using bare parameters of the Hamiltonian we conclude in agreement with the the conven- 
tional wisdom that disorder or periodic potential does not allow conserving current states 
in ID. The crucial point for the study presented here is that the SC in the interacting ID 
system are always "dressed" by density fluctuations. Any attempt to change / would result 
in driving the wavefunction of the whole ring into orthogonal state (that is (/ | I' ^ I) — > 
with macroscopic accuracy - see also below). Now, the renormalized amplitude connecting 
different numbers I is size-dependent for the lowest excited states. Depending on the in- 
teraction parameters between the particles in the original model, one can obtain arbitrarily 
long relaxation times for the topological excitations. 

In this paper we discuss the stability and relaxation of SC in ID interacting systems. To 
be definite, most of the time we consider Bose systems. As is well known ||, the statistics 
and interaction effects are indistinguishable in ID, so that fermions are exactly described 
as hard-core bosons, and spin chains are converted into the interacting bosons by Holstein- 
Primakoff transformation. Our main interest is in the low-energy excitation spectrum of the 
ring in the superfluid phase, including the energy splittings between the resonating SC and 
their hybridization with the sound modes. 

Experimentally SC may be observed in the response of the ring to the vector potential. 
Suppose that at the initial moment of time the system is in the ground state with some value 
of / = Iq- In the general case Iq is a function of the magnetic flux through the ring and 
the parity of the total number of particles in the system: for bosons and odd parity fermion 
systems in zero flux Iq = 0; fermion ring with the even number of particles is equivalent to 
a bosonic system with the half-quantum of the magnetic flux through the ring. In this case 
the ground state is degenerate between 1 = 1 and 1 = 0. By switching on instantaneously 
an integer number of flux quanta, <fi = rrKpo = mhc/e, we excite nonadiabatically the super- 
current state with I = Iq + m. One may then study the time evolution of I by measuring 
e.g. the oscillations and decay of the magnetic moment generated by the current. In a static 
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experimental setup the relevant quantity is the magnitude and sign of the persistent current 
near the points of the ground state degeneracy. In bosonic and odd parity fermionic systems 
these points are given by <p ~ (m + l/2)0o- Near these points the persistent current has 
the maximum possible amplitude and abruptly changes its sign. The width of the crossover 
region is very small | — (m + 1/2) O |^ <Po an d directly related to the frequency of the 
current oscillations in the dynamic experiment. In the last section we will describe in some 
detail the most promising experimental ID system - edge states of the 2D electron gas in the 
regime of fractional quantum Hall effect (FQHE). Here the SC are related to the electrical 
charge transfer between the edges in the Carbino disk geometry thus making it possible to 
study the phenomenon by monitoring the charging effects between the edges. Finally, we 
will demonstrate that the physics of SC relaxation in ID systems can be directly applied to 
a 3-dimensional superconducting ring with the ring length being much longer than the wire 
cross-section diameter. 

The other way to excite the SC is by pumping the system with the high-frequency 
electromagnetic radiation in resonance with the energy of, say, the first current state. 

The paper is organized as follows. In section [TT] we formulate the effective-action approach 
to the problem and show it to be very convenient for unified description of the critical 
points and renormalization group (RG) equations for the superfluid-insulator transitions 
in commensurate and disordered systems. Next we calculate the energy splitting, Aab, 
between degenerate current groundstates (in the presence of the half-quantum of magnetic 
flux or its equivalent) using an instanton technique. This splitting and its dependence on 
the ring length L is the main mesoscopic characteristic of SC. It gives the frequency of the 
coherent current oscillations in the superfluid phase, while its divergence signals about the 
superfluid-insulator transition. All calculations in section [TJ are done with the logarithmic 
accuracy. 

In section [ITT] we formulate an alternative approach based on the effective Hamiltonian 
for the SC dynamics. This way we calculate precisely the hybridization coupling between 
the topological states and phonons, which allows to fix all the ratios between the lowest 
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energy levels. The other goal is the possibility to account for the resonances between the SC 
and sound waves and describe correctly level crossings. All the results obtained analytically 
in sections [Tl] and [TIT] are tested in section [TV] by intensive numeric calculations done on 
finite-size Hubbard-like ring models. We present the low-energy spectra for various systems 
and discuss the procedure of their extrapolation on larger system size by combining the 
numeric data with the analytic RG equations. In section |V| we discuss possible experimental 
setups for the study of SC. We also considered in some detail two important experimental 
systems, i.e., the 2D electrons in the FQHE state and a long, narrow 3D superconducting 
ring, as the most promising systems where our results may be tested. 

II. EFFECTIVE- ACTION APPROACH 

Since the stability of the SC states is associated with the conservation of the topological 
invariant / ([!]) the most natural language for the description of its violation or even complete 
disappearance under various conditions is that which can directly take into account processes 
changing the value of /. In this section we introduce such a language in terms of the effective 
action for the phase field in (1+1) dimensions. 

Consider Popov's hydrodynamical action |3|] with the density-field fluctuations integrated 
out: 

S[*] = Jdxdr {m (x)$' T + U s m 2 + ■ (2) 

Here r) is the phase field (x is the spatial coordinate and r is the imaginary time), 
no(x) is the local number density with the coordinate dependence coming from this or that 
short-range external potential, A s is the superfluid stiffness, and k = (dfi/dn) -1 is the 
compressibility (/a is the chemical potential and n is the macroscopic density). As we shall 
see below, the dependence of no on x in (0) is of prime importance being responsible for a 
particular type of possible violation of superfluidity. So n (x) can not be replaced by n. 

Zero-point fluctuations with the change of I are represented by the topological defects 
(vortices) in the field r). Hence, the question of the existence in a ring system of the 
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metastable groundstate-like excitations with nonzero expectation values of / is the question 
of statistics of the vortices. If in the long-range limit the vortices are absent (bounded in 
short-range pairs) such states are possible. If separate vortices are present at the arbitrarily 
large scales one can not introduce the topological quantum number and the only possible 
equilibrium expectation value of / is zero. In this case the very description of the system in 
terms of the effective action (fj) is, generally speaking, no longer valid because of divergent 
long-range renormalizations. 

The self-consistent criteria for superfluidity under different conditions (commensurate 
external potential, disorder, or a special case with the single impurity) follow directly from 
the action (0). With rescaled variable r this action reads 

Sm = Jdxdy {m (x)$; + ^(V$ ) 2 } . (3) 

Here y = cr (c = JA s /k is the sound velocity), 

K - l = itJk~k . (4) 

The first term in ([3]) is of topological nature and "reacts" only on topological defects. Inte- 
grating over y and introducing 

rx 

7(2) = — 2n / n (x')dx' (5) 
Jo 

one may rewrite S[<&] as 

sm = ^{v<$>) 2 +iY.pn^), (6) 

where j enumerates the vortices, pj = ±1, ±2, . . . and Xj are the charge and the x-coordinate 
of the vortex j, respectively. The effective action (^|) has the form of XF-model (the 
first term) with the additional feature that each vortex brings in a spatially- and charge- 
dependent phase (the second term). This phase plays an extremely important role. To 
realize this, consider a homogeneous system. In the absence of the second term in the 
action (f|) the system would be equivalent to the XF-model and thus would demonstrate 
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Kosterlitz-Thouless transition at K — 1/2 which means destruction of superfluidity. How- 
ever, it is clear that the homogeneous system is always trivially "superfluid" because of the 
momentum conservation. Hence, no phase transition actually should occur, and it is the 
second term in Eq.(|6]) that corrects the result. Indeed, in the homogeneous case j(xj) is a 
linear function of Xj and the integration over Xj makes the net contribution of the vortex j to 
the partition function to be zero. A similar picture takes place in a periodic incommensurate 
potential (and, in particular, in a lattice at an incommensurate filling). In this case 7(x) 
may be represented as 



where x = si + £, / is the period of the potential, s is the corresponding number of periods, 
£ G [0, 1] is the reduced coordinate, v = Jq n (x)dx = nl is the filling factor. The second term 
in the r.h.s. of Eq.([FD leads only to a certain short-range renormalization of the statistical 
weight of the vortex and, thus, does not play an essential role, so we omit this term from 
now on. (Note also, that in a lattice this term is absent from the very beginning.) The first 
term in the r.h.s. of Eq.([?p is linear in s. This means that just as in the homogeneous case 
the net contribution of vortices to the partition function is zero, provided v is not a rational 
number. Hence, no phase transition occurs. 

In the case of commensurate filling the discreteness of s becomes important. Now there 
appears a class of vortices whose net contribution is nonzero. These are the vortices whose 
charges are the multiples of p, the denominator of the filling factor. All these vortices enter 
the partition function with one and the same spatially independent phase, so one deals with 
the XY-mode\ with the constraint that the charge of a vortex be a multiple of p. Thus, the 
system should experience a Kosterlitz-Thouless-type phase transition at K = p 2 /2, which 
can be identified with the known superfluid - Mott insulator transition in a commensurate 
system Note, that the critical value of K follows immediately from the effective 

action (||) by the free-energy-sign argument for a single vortex |5]]. This argument is exact 
and need not a resorting to the RG, provided K is understood as the long-range asymptotic 

7 




(7) 



quantity |]]. Nevertheless, RG equations are of interest for some applications (see below). 
They read (cf. 0) 

dK/dX = w 2 

(8) 

dw/dX = (2 - p 2 /K) w 

Here A = InR, R = y/x' 2 + y' 2 being the characteristic scale of distance, K(X) is the meso- 
scopic value of K. Eqs.(§) are the standard Kosterlitz-Thouless RG equations written 
for the vortices of the charge p. As usual, (w(X)) 2 is proportional to the number of vor- 
tex pairs of the size ~ R in the area ~ R 2 and is much less than unity in the region of 
applicability of Eqs.(^j). Mesoscopic values of A s and k are related toK(X) by 

c being irrenormalizable quantity. 

Considering disordered superffuid system it is convenient to rewrite the second term in 
(|) in terms of the pairs of vortices: J2jPjl( x j) ~> X< 1c, 7( x iC' x 2c)> where ( indexes the 
vortex pairs, = 1, 2, 3, . . . is the charge of the vortices in the pair (, and a^c are the 
x-coordinates of the positive and the negative vortex in the pair (, respectively, and 

f X2 

7(^1, X2) = — 2n / n (x)dx . (10) 

J XI 

Being interested in the long-range contribution of a vortex pair we may average its phase 
factor over short-range translations along the x-axis: 



e 



i2irq ; y(x 1 ,X2) _^ ( e i27r< ? 7(^1+^.^2+0 \ _ f ( Xl — £ 2 ) (11) 



Since no{x) in a disordered system is a random function with some microscopic correlation 
length, the function f{x\ — x 2 ) decays at microscopically small | x% — X2 | and may be 
replaced by oc 5(xi — x 2 ). We see thus, that in a disordered system the second term in 
the effective action (^|) is equivalent to the constraint that only "vertical" vortex pairs 
be present. This constraint changes the entropy per vortex as compared to the regular 
commensurate case. Now the entropy per vortex is (3/2) In L+ finite term since the entropy 
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of a pair is obviously 3 In L + finite term. The energy of a vortex with the minimal charge is 
K~ l In L + finite term, so the free-energy-sign argument yields the critical macroscopic value 
of K = 2/3, in agreement with the result of Giamarchi and Schulz || obtained from other 
considerations. The renormalization group for the disordered case may be easily derived 
from the effective action in analogy with the Kosterlitz-Thouless treatment ||. The new 
feature associated with the above-mentioned constraint is that renormalizable quantities 
now are A s and c while k can not be renormalized by the vertical vortex pairs. The RG 
equations read now 

| dK/dX = Kw 2 

(12) 

I dw/d\= (3/2-l/K)w 
Here A and w have the same meaning as in Eqs.(§). Mesoscopic values of A s and c are 
related to K(X) by 

A * (A > = ^¥(W < x> = Jm- (13) 

Since in the disordered case there takes place a renormalization of c, a question may arise of 
what quantity c should be used for the introduction of the variable y = cr. For the derivation 
of Eqs.flT2|) it is natural to choose c to be some small-scale value of c(A). For the application 
of RG to a particular mesoscopic system the most appropriate choice is c = c(A*), where A* 
is the large-scale cutoff for A, related to the size of the system. We notice, however, that 
this difference leads only to small corrections in the parameter w 2 <C 1 to the RG equations, 
and thus may be neglected. 

Finally, consider a special case of a regular system with a single impurity (or any other 
defect). As it follows from the above analysis for regular and disordered systems, in this 
case the effective action corresponds to XY-mode\ with the constraint that the vortex pairs 
be vertical and located only in the vicinity of the line x = Xq, Xq being the coordinate of the 
defect. The free-energy-sign argument (entropy per vortex is In L + finite term) immediately 
yields the critical value of K = 1, in accordance with the result of Kane and Fisher [Q. 
In the case of one defect no renormalization of K occurs. Physically it is clear since the 



transition is associated only with the renormalization of the strength of the defect: effective 
tunneling constant tends to zero at K > 1 ||. 

Now we use the effective action to consider the splitting Ai of the SC level / in a finite- 
size system as a function of L and K. This splitting is due to the tunneling between the 
states | I ) and | —I ) which means that the true eigenstates of the finite-size system are 
superpositions of 1 1 ) and | —I ) . To calculate Aj ab initio one may take advantage of the 
conventional instanton approach starting from the formula 

( -/ | e iHt | I) = e - iiE ° +El)t sin(A J t) , (14) 

where E g is the groundstate energy, Ej = 2tt 2 I 2 A s /L is the supercurrent energy, H is 
the system Hamiltonian. The next step is to use the imaginary time t = —i(3 in order to 
represent the l.h.s. of flT4]) in terms of the effective action (0). Since the functional integral is 
naturally introduced in terms of the coherent states, we represent | / ) as the expansion over 
the set of coherent states {| ip )}. Actually, we are interested only in the relation between the 
expansion for 1 1) and that for the groundstate | ). In the homogeneous case this relation 
is trivial because of the Galilean invariance: If 

|0) = Jd^Q[^]\^), (15) 

then 

\I) = |£tyQM |e <2 ^V>- (16) 

In the general case one may think of the state 1 1 ) as being obtained from the homogeneous 
one (|i~6|) by semiadiabatic turning on the external potential. The word " semiadiabatic" 
means that the turning-on time is much greater than the inverse interlevel separation oc L, 
but is much less than Aj 1 . Integrating over all the variables except for the long-range phase 
field $ we have 

( - J | e~ Hp \I) = A J e~ 5[#1 , (17) 
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where A is some independent of L and / normalization constant, the field $(x,r) is deter- 
mined at < r < P and satisfies the boundary conditions 

[ L dx -^-§(x,0) = 2ixl , f L dx ^-$(x 7 (3) = -2nl . (18) 

•J CJJu J CJJb 

To the first approximation, semiadiabatic turning on the external potential does not affect 
Eq . (|l7l) , provided (3 is much greater than the turning-on time. 

Calculation of the integral in (|T7|) may be performed by a shift of the variable $: 

$ = $i + $o, (19) 

where the shifting field $o is chosen to satisfy (below we turn to the variable y = ct and 
correspondingly = c/3 ) 

A$ = (20) 

with the boundary conditions ([Lj) ($ — > $ ). Then £>[<&] = 5[$i] + S[<3> ] and $ a may be 
integrated out taking into account that 

A J D<$> 1 = ( | e~ Hp | ) = e"^ . (21) 

As is easily seen, to meet the boundary conditions $o should contain a number of vortices 
(=instantons) with the total charge 21. So besides the integration over $ x there should be 
also an integration over the positions of the instantons: — > D&i Hj^dxjdyj, (xj,yj) 
being the position of the vortex j. Generally speaking, there is also summation over the 
number of instantons, iV*. For our purposes, however, we may confine ourselves to the case 
(3 <C Aj 1 , where the main contribution comes from the N* corresponding to the minimal 
possible value of £>[$()] ■ This implies that all vortices are of the minimal possible charge of 
the same sign. We notice also that in this case sin(— iAj/3) ~ —iAj/3. Consequently, from 



and (pip we obtain the relation for A/: 

n 

3=1 
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J J] dx 3 dy je - s[ * o] oc e- ElP Aii3 , A//3 < 1 . (22) 



Consider the most simple case of the splitting of the two-fold degenerate groundstate 
level in the system with the gauge phase 7T, which we will refer to as Aharonov-Bohm (AB) 
splitting, Aab- Upon the elimination of the vector potential by the gauge transformation 
this case corresponds to / = ±1/2. It is worthmentioning certain ring systems in which such 
a case occurs automatically: fermions with the even number of particles; lattice spins with 
antiferromagnetic interaction and lattice bosons with the positive hopping amplitude, both 
at the odd number of sites. 

For 1=1/2 minimal S[§o] corresponds to $o with the only one vortex of the charge 
equal to unity: 

S[%] = * 7 (x ) + + K~ l In L + o(L _1 ) , (23) 

where x$ is the x-coordinate of the vortex core. Taking into account that 

const one defect 
~ \f~L disordered system 



o 



L dx e n{xo) oc < 



(24) 

L commensurate with integer v 



otherwise 

from (|22~|) we obtain the estimate for the relative value of the AB splitting. 

1/ L i/k-i one defect 
l/L/VK- 3 / 2 disordered system 



oc < 



E 



AB 



(25) 

\jL x l K 2 commensurate with integer v 



otherwise 

Here Eab = -E/=i/2- Note that in the disordered case the value of splitting is sensitive to the 
particular realization of disorder and experiences relative fluctuations of order unity from 
system to system. Also, in the regular incommensurate case as well as in the commensurate 
case with a fractional filling factor Aab is equal to zero because of the exact cancellation of 
contributions from instantons with different Xo's. 

For disordered and commensurate systems the laws (|2"oD, generally speaking, take place 
only at sufficiently large L, where K reaches its macroscopic limit. So they need be revised 

12 



in the vicinity of the critical points, where renormalization of K is essential and, moreover, 
is very slow with respect to L. (Such a problem does not emerge in the single-impurity case, 
since there is no renormalization of K.) To allow for the renormalization we notice that the 



only place in the effective action ( |23"D where smaller scales of distance may be important is 
the vortex core energy Ey(L) = K~ 1 \nL, which now is to be estimated more accurately. 
This may be done just by noting that according to the above-mentioned meaning of the 
variable w in RG equations it is related to Ey by 

( L 3 e ~2E v (L) disordered 
(26) 
£4 g -2£V(L) commensurate 

Thus we obtain 

^ « w{L) , (27) 

where w(L) should be taken from RG equations. In particular, at the critical point, where 
both (D and © yield w(L) ~ 1/lnL, we have 

— — oc 1/lnL (at the critical point) . (28) 
Eab 

In the case of commensurate system with a fractional u oc 1/p only the instantons whose 
charges are multiples of p are relevant. Clearly, this leads to the fact that only some special 
SC levels are split. Namely, the levels with J's being the multiples of p/2. Consequently, at 
an odd p split SC levels arise only at the gauge phase tt, while at an even p the gauge phase 
7r renders all SC levels degenerate. Estimation of the splitting of the first split SC level 
(/ = ±p/2) is done in the complete analogy to the estimation of Aab- The only difference 
is that the charge of the instanton now is p rather than unity. So we readily get 

— — oc 1/L p2 / K ~ 2 commensurate with v oc 1/p . (29) 

Ep/2 



This result is a generalization of Eq . (|25|) for the commensurate case. Obviously, Eqs. (p7"U28D 
(with AB — > p/2) are also valid. 
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If we try to estimate by the same technique the splittings of higher SC levels, when more 
than one instanton is necessary, we have to deal in ( |22| ) with a divergency of the integration 
with respect to y/s. This is a standard problem one faces when considers in imaginary time 
the tunneling between the excited states rather than between ground-state ones. Instead of 
thinking of a trick of circumventing this difficulty, we prefer to turn to another language of 
description. In the next section we introduce the effective Hamiltonian, whose relation to 
the effective action (§) is clearly traceable, and which allows to treat the problem of the SC 
splitting as well as the hybridization with phonons most naturally. 



III. EFFECTIVE-HAMILTONIAN APPROACH 

In this section we demonstrate how one may construct an effective Hamiltonian for the 
SC dynamics and calculate the ratios Aj/Aab as well as the resonances with the phonons. 
In the Fourier-transformed fields the Hamiltonian of the superfluid phase can be written 

as 

H = Y,u q blb q +^-(v N 5N 2 + 4v I I 2 ) , (30) 

q ZL 

where vi = 7rA s , vn = (kk)* 1 , b\ is the standard Bose creation operator, u q ~ c | q | in 
the long- wavelength limit. For completeness, we have introduced the term associated with 
the energy change when we add particles to the system, 5N — N — No- 

We start by noting that the instanton is nothing else but an effective vertex connecting 
nearest SC. Thus in the most general form we may present the instanton at the point x as 

Hint{x) = £ (^e***) 0(b) c\ +1 cj + /i.e.) , (31) 

Here c\ +1 is the raising operator of the quantum number /, 7(0;) is given by Eq.(f|), and O is 
an operator acting on the phonon states. The dependence on the coordinate comes from the 
fact that an intrinsic momentum of the current state in the homogeneous system equals to 
I2irn and the Hamiltonian (|31|) changes / by 1. Indeed, the SC is obtained by shifting all the 
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particles in the momentum space by / quanta with the total momentum change NI 2n/L. 
In disordered systems the particle density is no longer homogeneous and we must replace 
the average particle density with the local one n — > no(x) according to (|5|). To satisfy the 
RG equation for the effective vertex w (see also below) the scaling dimension of O should 
equal to 1 — 1/K. 

As shown by Haldane |§ (and will be demonstrated explicitly below for the particular 
case of electron backscattering from the impurity) the bosonic field responsible for the current 
relaxation is 



6{x) = (0„ + 7(*)) +iiE J^ 7 sign(g)(6 (? + bl q )e~^ . (32) 

v K g^ V I Q I L 

and the original particle field is given in terms of exponential functions of 9{x). By noting 
that the phase 8q is the canonical momentum for /, the resulting effective Hamiltonian may 
be rewritten as 

H mt (x) = £ C^e^-HW) 0(b) 4 Q c eo + h.c) , (33) 

-K<0 O <-K V L J 

Clearly, the first two terms in the exponent are those in the definition of 9(x) given above, 
and we identify the operator part of the vertex O in Eq.([H]) with the q ^ components of 
the phase field 9(x) 

6(b) = -enviiB^x)} , (34) 

where the first factor is necessary for the correct scaling of O. The sum over momentum 
has an upper cutoff at some scale q < 2tt/A where one has to define the "bare" amplitude 
Va which depends on the microscopic parameters of the initial Hamiltonian and the cutoff. 
One can not obtain the exact value of the parameter Va from this analysis (except for the 
scaling behaviour mentioned above). The other way of deriving the effective Hamiltonian 
mathematically rigorously is by mapping the effective Lagrangian Eqs. (0fj) to the equivalent 
sine-Gordon model 0, which is in fact done by Eq.(|33"D since Hmt ~ cos0(x). 

There is yet another simple microscopic derivation of Eq.fl54j) for the particular case 
of electron backscattering from the impurity potential. In the bosonic representation the 
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original right- and left-moving electron fields (labeled by the index p = ±1) have the form 



¥ p (x) = _L Up e-^exp ^^sign(g) C pq (b q + bl q )e 



■iqx 



(35) 



^pq 



cosh(/\~) , pq > 
sinh(/\~) , pq < 

where the Fermi momentum fcp equals to ttN/L, and U p is the ladder operator increasing 
the number of p-particles by one. If we substitute these expressions into the backscattering 
Hamiltonian, Hb = Vb^I^i + h.c, which is responsible for the current relaxation on im- 
purities, we immediately arrive at the effective Hamiltonian Eq . fl33| , |3"4D [|13|1 . An important 
new point here is that it is valid for any ID system in the superffuid phase, provided one 
uses the proper definition of the quantum number J, which has the particular meaning of 
the difference between the number of left- and right-moving particles for fermions only 
A short physical discussion is in order here. One may consider the effective Hamiltonian 



(|3T|) as resulting from the coupling between the SC states and sound waves of the form 



Ic 



K q^O V L 

As is well known, this interaction can be eliminated by the standard on-site polaronic trans- 
formation in expense of the additional operator structure O in the hopping term. Thus 
any SC state causes a long-range deformation in the phonon field, and this deformation is 
growing linearly with the quantum number I. In the fermion system this is equivalent to 
"dressing" each extra left/right-moving particle by density fluctuations. Now any attempt 
to change the SC state will result in the readjustment of the phonon vacuum to the new 
value of /. It turns out that the dominant contribution to the effective low-energy vertex is 
due to this reconstruction of density fluctuations rather then due to the "bare" process of 
current relaxation which defines the high-energy prefactor here. As K goes to zero the cou- 
pling between the SC state and density fluctuations diverges. This limit describes the case 
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of weakly interacting bosons and quasi ID superconducting rings (see section |V|), where even 
the lowest SC state (which still has a finite momentum!) involves large number of particles 
with small momenta moving in the same direction. Now the process of current relaxation 
is through the backscattering of many particles simultaneously, rendering the low-energy 
vertex extremely small. Thus a one-particle description of SC in terms of the number of 
left- and right-moving carriers is clearly misleading for K — > 0; the relevant language in all 
cases being that based on the quantum number I. 

One may now straightforwardly check that the renormalization group equations in the 
disordered case are correctly described by the effective Hamiltonian. The RG equation 
for the vertex is obtained by a simple averaging of the operator O over the momentum 
range 2it/A' <\ q \< 2ti/A where A' = A + dA. In the second order perturbation theory 
the renormalization of the superfluid energy is due to the virtual transition of the system 
to the neighbouring current state with the simultaneous emission of the sound mode with 
27r/|g|e(A,A') 

1 f dq 



5E I = -{V/Lf-2 f 

K J(A,A' 



(A,A') q 



1 

+ 



(37) 



,2-kc/A + E I+ i — Ei 27rc/A + - E T _ 
where Ej = 2nviI 2 /L. By expanding the denominators up to the second order in 1/L we 
find the RG equation for the current velocity: 

dlnvi W 2 A 

din A % 2 c 2 K 2 L ' 1 ' 

The renormalization of the sound velocity is due to the virtual transitions from I to / ± 1 
with the simultaneous absorption of the sound wave and emission of the high-frequency 
mode. In full analogy with the previous calculation we have 

dine 2V 2 A 

din A n 2 c 2 K 2 L ' 1 ' 

As expected, one may neglect the renormalization of Vj and c in the superfluid phase when 
dealing with the single impurity scattering (the r.h.s. of Eqs.( |38| , |3"9|) is 1/L in this case). 

For disordered potential the effective Hamiltonian is an obvious sum over the impurity 
positions, and in Eq. (|3~l~D one has to replace V Y, Xi W 7 ^ ~ V \Lj 'A) 1 / 2 . Now the 
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scaling equation for the vertex is changed to dV/d\nA = V(3/2 — 1/K), in agreement with 
Eq.(O). Correspondingly, in Eqs.(|38|,|3"9|) one has to replace V 2 by V 2 L/A. By introducing 



the dimensionless coupling parameter w = 2V/ncK we rewrite the RG equations as 



d In Vj 2 
- — w 



din A 

or, using the universal relations c = vjK = v^/K 

dK 



din A 

d In c 9 . , „ 

-w 2 /2 , (40) 



din A 

Thus we have reproduced the RG equation (|l|) of section ||. 



w K ; Vjsf — const . (41) 



In exactly the same manner one can construct an effective Hamiltonian for the commen- 
surate system. Let n — I /p. Then the vertex connecting the states I and I ±p is defined 
by 

dx (V 



H coram 



? / f {j^ 6 ^ c Uci + h.c) . (42) 



Due to the translational symmetry of the commensurate problem there is no renormalization 
of the sound velocity For the same reason only states with zero total momentum are 
hybridized with the topological excitations and e.g. the correction to the superfluid energy is 
defined by virtual transitions to the state with two high-frequency modes q\ = —q 2 G (A, A'). 
Now 

dK 



p /16u> ; c = const . (43) 



dlnA 

The coefficient in front of the w 2 term can be absorbed into the redefinition of w (cf. Eq.([|)). 

The effective-Hamiltonian approach reduces the discussion of the superfluid - insulator 
transition to that of a "particle" localization at a given site /. In the ground state (when 
A = L) the current Hamiltonian is identical to the tight-binding model with the hopping 
amplitude Vl/L in the static external potential given by the supercurrent energy 

H{L) « E - 0M,) 2 c\ Cl + c ^ + h - c ) ■ ^ 



If under scaling one finds Vl <C fj, then the potential energy prevents the current state 
from derealization, and we have perfect metallic conductance. Only at some special values 
of the magnetic flux through the ring (e.g. integer or half-integer flux), when the states / 
and —I are in resonance, the current is delocalized between the two states, but the energy 
splitting of the resonance is much less than 2ttvj/L. 

The superfluid-insulator transition is described as a derealization of the current state, 
when Vl ^> Vj. In this regime the topological invariant is delocalized in a band and / is no 
longer a well-defined quantum number Jl3| . 



Let us now calculate the splittings between the resonating states. If we neglect the 



hybridization with phonons and work within the Hamiltonian (|44), then the splitting is 
given by the perturbation theory of 2/-th order in Vl 

y2I I-l I 

A 7 = 2^- TT 

*EL (lL-\ 21 1 
L \2nvj) [(2/-1)!] 2 ' 1 1 



If the ring is in the flux <fi = <fio/2 then the splitting between the states I and 1 — I is 
described by the same expression (]45f) with the replacement I — > I — 1/2. 

In deriving Eq.(^5|) we considered only the path from / to —I which corresponds to the 
ground state of the bosonic field in all intermediate states. In the general case one has to 
consider the possibility of exciting/absorbing an arbitrary number of sound waves in any 
intermediate state. To this end we slightly " undress" the ground-state vertex and write the 
final low-energy Hamiltonian in the form 

Heft = E 2 -^(i - 0M>) 2 c\ Cl + E (i^y c '+^ + h - c ) ■ ( 46 ) 

which is cutoff independent. First, let us consider the case of single impurity, when the 
spatial dependence of Vl is S(x — Xj) at the impurity position Xj. Using standard properties 
of bosonic operators we readily calculate all the matrix elements 

(TO I \{N q + S q }) = Ur(N q ,N q + S g , T ^ l ), 
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x S/2 oo (_l)n( N , Q , \\ 1 S > 

r(N,N + S,x)= / 6X en, ~ ■ ( 47 ) 

In principle, the problem of the low-energy spectrum is solved now, because we know ex- 
plicitly all the matrix elements. In a large system the coupling Vl is macroscopically small, 
and apart from extremely narrow resonances at 1/K = integer, when one has to solve a 
finite-size secular equation, it is sufficient to restrict ourselves to the perturbation theory in 
Vl. Below we calculate analytically the energy splitting of the first SC state. 

First, we note that the coupling to phonons does not modify the splitting of the ground 
state in a ring with = 0o/2, simply because there are no intermediate states in this case 
(see Eq.@)) 

4ttvi 



iV L /2nvi) . (48) 



L 

With the coupling to phonons taken into account the expression for the first SC splitting 
acquires the form 



L _ \r(0,S q 



2n 



Next we introduce the integer variable M and write down the spectral representation for 

the 5-symbol 5(M — J2 q (L \ q \ /2n)S q ). This allows to calculate the contribution of each 

phonon mode separately 

_ 47TO/Ai r* dt ~ e~ iMt { 2 ~ e ipt 

~ LA\ B ~ Uto ^ Q 1-KM e ^\K^[V 

oo df e -iM t{l _ ^ty/K 

~~ £±oJ-*2ic 1-KM ' ^ ' 

Straightforward evaluation of this expression gives the final answer 

R - «M m + 2/K ) 

sm(7r/K) T 2 {l + 1/K) ' 1 ; 

Note that R = 1 if we neglect the hybridization with phonons. 

Thus we found that for 1/K — integer the splitting diverges due to the resonance between 
the SC and phonons. For small Vl this divergence is accurately described by Eqs. ([50"1 , |51~1) 
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except for a tiny region of order Vl/vi around integer values of 1/K. However in this region 
the ratio R is already very large, and the dominant contribution comes from the coupling 
to a small group of phonon states with the same energy ss 2irvi/L. For the particular 
case of 1/K ss 2, which will be compared in the next section with the exact numerical 
diagonalization, one has to take into account the states | N q — 2), | N_ q = 2), | N q — 
l,N- q = 1), | N 2q = 1), and | N_ 2q = 2) where q = 2n/L. The hybridization matrix 
elements between these states and SC squared are 2, 2, 4, 1, and 1, correspondingly. By 
writing a secular equation for this group of levels we obtain a complete picture of level 
crossings between the SC and sound modes. 

Quite generally at the point of resonance all energy splittings are of the same order 
A/ ~ A AB . Indeed, for integer 1/K the instanton is always between the states with the 
same energy (for the combined system SC + phonons), so the splitting is given by the 
amplitude of a single transition Aab no matter how many instantons are required to connect 
the resonating SC states. 

This technique may be applied for calculating the energy splittings in commensurate 
systems as well. An important difference with the impurity case is that now the supercurrent 
energy change is SEj_ p j = 2irvj(2pl —p 2 )/ L and resonating states are given by / = p/2±mp, 
and only sound waves with the zero total momentum are hybridized with the SC. We find 
that for commensurate systems we always have a resonance between the SC and sound 
at the point of superfluid-insulator transition. Indeed, for K = p 2 /2 the energy change 
5Ej± p j = 2nc/L 2(2n — 1) is given by even number of the phonon energy quanta. 

Strictly speaking, there is only one cutoff-dependent parameter Vl (apart from the Lut- 
tinger Liquid parameters c and K) which controls the low-energy spectrum of the system. 
Once this parameter is fixed, say from the Aab splitting, all higher energy splittings and 
level crossings can be found explicitly However in the disordered system we face a problem. 
By construction the low-energy vertex Vl(x) in this case is a random function of coordinate, 
and depends on the particular disorder realization in the sample. Now, even if we define the 
mean value of < Vl > from the lowest splitting Aab, we still are uncertain about the inte- 
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grals / L dxe ±l2lvmx I l Vl(x) which define the hybridization matrix elements with the phonon 
modes and take care about momentum non-conservation in disordered system. The only 
thing which is certain is that for the lowest levels with m ~ 0(1) these integrals have random 
values of order < Vl > in magnitude. It makes the direct comparison between the theoretical 
calculation (|5l|) and the numerical result for R inconclusive, unless additional parameters 
are introduced into the theory to account for the random variations of the above-mentioned 
integrals (see more in the discussion of the numerical results in section 



The above discussion dealt with the properties of an isolated ID ring. Obviously, by 
considering the lowest energy levels we lack any kind of dissipation in the system. In the 
ideal picture the current state I obtained, e.g., by nonadiabatic switching on the magnetic 
flux through the ring will persist forever, or demonstrate undamped oscillations with the 
frequency Aj. The most effective dissipation mechanism is due to the coupling between 
the ring and substrate. One has to take an insulating substrate in order to eliminate the 
conduction electrons which form a thermal bath with the largest low-frequency density of 
states, and thus are the most dangerous source of decoherence. To couple the SC states 
to phonons of the substrate we have to consider the variation of H e ff originating from 
the lattice long wave-length distortion u(x). The effect of u(x) is two-fold: to change the 
kinetic and potential energy. For example, the distortion with the wave-length of order L 
(which can absorb the energy difference 5Eij_ p ) may influence the total length of the ring 
5L ~ Lu and result in the u-dependent terms both in the kinetic and potential energy. 
For disordered sample the positions of impurities {xj} depend on the substrate deformation 
giving SVl ~ Vl,u. One may also take into account the dependence of K (or the particle- 
particle interaction) on the inter-particle separation. 

It is easy to check that both the nondiagonal and diagonal couplings give the same value 
of the relaxation time r. We write the diagonal interaction Hamiltonian as 

Vi Vl 



H diss ~ o:—u ~ a- 



(52) 



Y,\lk/Ms{d k + d\) 
. k 

where d\ are the bulk phonon creation operators, a is the dimensionless parameter depending 
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on the microscopic details and geometry of the particular system, M is the substrate ion 
mass, and s is the bulk sound velocity Now, using standard results for the relaxation time 



in a double- well system coupled to phonons (see e.g. we find 



-'-(W^-^. ( ^ 0) - (53) 

Here p is the substrate density. If we couple the substrate to the nondiagonal terms in 
Eq.(f46l) by expanding 

V L /L -> V L /L (1 + a'u) , 

then the final result for r _1 coincides with the Eq(|53D. 

The relaxation rate should be much less than the level splitting Aj if we want to observe 
the coherent oscillations. In the mesoscopic limit we find Ajt ~ V 2 ^L 4 . Substituting 
here the scaling behaviour of Vl, we obtain the criteria defining which levels will demonstrate 
coherence at T = 



/ < < 



(1 + K)/(l — K) single impurity 

(2 + K)/(2-3K) disordered • (54) 

P 2 /(p 2 ~ 2K) commensurate 

We see that the broadening of the two lowest doublets Aab and A x is always small. Very 
close to the superfluid-insulator transition we can safely neglect dissipation for all levels. To 
get a better feeling of how long the relaxation time could be in a real system, we estimate 



t 1 from Eq.(53) using 7~1, Vi ~ s, L~ 10 4 A°, and the maximum possible Vl ~ v\. For 
this particular choice of parameters r ~ pL 5 ~ 1 week. 

IV. NUMERICAL SPECTROSCOPY 

In this section we report exact numerical spectroscopic results for ID finite-size boson 
Hubbard models. The Hamiltonians we examine can be written in the most general form as 
(all the parameters are scaled in the units of the hopping amplitude) 
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f U 1 

# = \ -( a t a i+i + n - c + e * n i + iz n i( n i - 1) + f • (55) 

i=i ^ ^ J 

Here aj is the annihilation operator on the site i\ rij = a+aj, and {ej} are random numbers 

uniformly distributed over the range from — W/2 to W/2. (In the single- impurity case all 

the 6j's are equal to zero except for ei.) The ring topology implies that the sum of the site 

indices is understood modulo N a , the total number of sites. The gauge phase ir, if necessary, 

is introduced by changing the sign of one of the hopping amplitudes (anti-periodic boundary 

condition) . 

For the calculation of the scalings of SC splittings with increasing N a we use hard-core 
models (U = oo) which at a given N a have much less Hilbert-space dimension as compared 
to the full models (U ^ oo), and thus allow consideration of larger systems. In this case the 
term with V is necessary to control the parameter K. In the case of a full model, where 
this term does not introduce any qualitative difference with respect to the on-site term U, 
we set V — 0. 



Our diagonalization procedure |22| is arranged as follows. First, starting from some trial 
wavefunction we calculate a certain number of the lowest energy levels and corresponding 
wavefunctions by modified Lanczos method [|l^|. The set of approximate eigenfunctions is 



reconstructed from Relay's tridiagonal matrix |[20|| , and the trial wavefunction is expanded 
in it. As is known, the set inevitably involves a substantial number of spurious states, 
due to numerical errors. These states, however, may be easily identified by their negligible 
contribution to the expansion of the trial wavefunction. Upon exclusion of the spurious states 
the set is subjected to the orthogonalization and correction by Newton method. The relative 
(with respect to a characteristic interlevel spacing) errors of the energy level calculation are 
typically of order 10~ 13 -v- 10~ n for the groundstate, and of order 10~ 9 -j- 10~ 5 for some ten 
first excited states. The problem of degeneracy is solved by repeating the procedure with 
a new trial function, chosen to be orthogonal to all previously obtained eigenstates. The 
combined set of eigenfunctions is again orthogonalized and corrected by Newton method. 
In the case of the regular system we take advantage of the translational symmetry to 
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proceed separately for each momentum sector. 

Now we turn to the results. To get an idea of how the spectra look like in this or 
that case, we show the pictures of spectrum evolution with increasing on-site disorder (Fig. 
1) and on-site interaction (at a fixed disorder) (Fig. 2), and as a function of the nearest- 
neighbor interaction in the case of commensurate hard-core model (Fig. 3). All the energies 
are reckoned from the ground states. Besides, in the case of disordered system we found it 
convenient to normalize the energy values in such a way that the middle of the first phonon 
doublet always corresponds to unity. With this normalization the value of the first SC 
coincides with K~ l . The level identification is performed in accordance with [p2 |, where the 
numerical spectra for a disordered ID ring were reported for the first time. 

Though it is not directly related to the subject of the present paper, it is interesting 
to discuss the behavior of the hard-core model at large negative Vs. From general con- 
siderations it is clear that at a certain negative V = V* there should occur an ordinary 
condensation into a macroscopic drop. So that at V < V* the first N a — 1 energy levels 
(associated with the motion of the drop as a whole with momenta = 1, 2, 3, ... , (Na — 1)) 
should lie very close to the ground-state one (corresponding to the zero-momentum motion 
of the drop). It is precisely what is seen in Fig.3. From the symmetry of the hard-core model 
(isomorphism with the spin- 1/2 system) it immediately follows that in our case V* = — 2 
is the point where, speaking the spin language (isotropic ferromagnetic Heisenberg Hamil- 
tonian), the transition to the ferromagnetic ground state takes place. To minimize the 
energy below this point a system with fixed ^-projection of the full spin (fixed number of 
particles in our case) would desire to be divided into two parts: with spins up and spins 
down, respectively (a part with empty sites and a part with completely occupied ones). 
In a regular system this symmetry breaking is hidden by translation invariance, but may 
be revealed by calculating some relevant correlator. Consider e.g. Q 3 = (n^n^n^) / (n^) 3 , 
where | i\ — 12 \ = \ i-2 — is | = | 23 — h \— N a /3. In the macroscopic limit we have Q3 = 1 for 
homogeneous liquid and Q3 — > for a condensed one. Our calculations for N a = 18 with 
the filling factor v = 1/2 yield Q 3 « 0.90 at V = -1.8 and Q 3 « 0.14 at V = -2.3. 
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In connection with the transition at V = —2 in the hard-core model it is worthmentioning 
that just at this point there exists an analytical Bethe-ansatz treatment recently reported 



by Sutherland pi| , which reveals the appearance of macroscopic bound complexes. 

The scaling of the single-impurity AB-splitting (Fig. 4) is in excellent agreement with 
the macroscopic theory. Here and below the dashed line ("theory") is plotted in such a 
way that it passes through the last (maximum N a ) point and has a slope corresponding 
to Eqs.(|25|, |2~9"D , the value of K being obtained from the numerical spectrum at this point 
(sound velocity can be taken from the first phonon level while A s is available e.g. from the 
first SC). 

The analogous results for the disordered system (Fig. 5) demonstrate a good agreement 
with the theory only at larger N a . In fact, this is what one might expect to get bearing in 
mind that the value of K we deal with (~ 0.5) is not very far from the critical one (= 2/3), 
so that the renormalization of K due to the disorder may be significant. In accordance 
with Eq.(^) the decrease of Aab with increasing iV a at K ~ 0.5 is rather slow. Hence, to 
reveal it against the statistical dispersion, characteristic of the disordered case, an extensive 
averaging is required. 

Fig. 6 shows the scaling of Ai in the commensurate case. The results at V = —0.5 
are in excellent agreement with Eq.(^), while at V = —1.5 there is a slight deviation 
(which, however, does not exceed 15%). Since this deviation has a tendency to decrease 
with increasing N a , it may be accounted for as a small-size effect. 

Fig. 7 shows the scaling of Ai for the commensurate system at the critical point and 
demonstrates the excellent agreement with the prediction of the macroscopic theory that 
the plot will approach a straight line as N a — > oo. It should be noted that for the half- 
filled ID hard-core model the critical point (corresponding to V — 2 where the model is 
isomorphic to the isotropic Heisenberg antiferromagnetic spin-1/2 chain) was extensively 
studied |[23| . In particular, it was shown that some excited levels demonstrate logarithmic 
(~ 1/ln N a ) relative shifts of their positions as compared to the macroscopic limit. Speaking 



superfluid language [p5|| , we understand that these are the states that have either non-zero 
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topological quantum number /, or extra particles/holes. (Purely phonon excitations can not 
demonstrate such shifts since the sound velocity is irrenormalizable quantity) 

To verify the correctness of the effective-Hamiltonian approach Eqs. (|3~T| , |3"4|) we calculate 
numerically the relation between the SC splittings Ai and Aab for the disordered (W = 0.25) 
system of 7 bosons on 11 sites. The coupling parameter U varied in the range 2.6 < U < 6.2 
where the first SC level crosses the group of phonon levels with energies ~ 4irc/L. The 
numerical ratio R = 2EiAi/Aab is shown in Fig. 8. Unfortunately, the direct comparison 
with the theoretical expression (0) is not very conclusive here because it was calculated 
by assuming that all 5 phonon levels with energies ~ Atcc/L are degenerate. However 
in a very small ring this is not the case even approximately as is clearly seen in Fig. 2. 
The single resonance is split and in the parameter range U < 3 other multi-phonon states 
(also strongly split) practically wash out any sign of the phonon spectrum quantization 
in multiples of 2nc/L. Under these conditions one may not expect to see any reasonable 
agreement between the calculated ratio Eq.(|5T|) and the numerical results on such small 
systems. 

One may try to improve the theory by adjusting the calculated ratio Q5TD to the actual 
positions of the phonon levels "by hand", e.g., by applying the procedure 
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(56) 



where the sum is over any group of levels characterized by {q, S q } with the matrix elements 
given by (0) (more precisely the ratio has to be calculated from the secular equation for 
the same group of levels). This looks like a reasonable approach to the single- impurity 
case. Still, in the disordered case we have further uncertainty in the values of the matrix 
elements because of the integrals / dxVi(x) exp{±ry(:r) ± i27rmx/L}, which depend on the 
particular disorder realization. Obviously, a single value of / dxVi(x) exp{i , y(x)} taken from 
the amplitude A^ is not sufficient to fix all these integrals for different phonon modes. In 
the general case we have to treat these integrals as random values, which for the lowest modes 
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are of order (Vl) in magnitude. The avoided level crossing in Fig. 2 tells unambiguously that 
all the symmetries are violated, and SC is coupled to the left- and right-moving phonons 
with different matrix elements. On the other hand we see clearly in Fig. 8 that the splitting 
Ai has extended resonance structure and the ratio R is typically much larger than unity in 
qualitative agreement with the theoretical prediction Q5T|). 

Finally, we touch upon a problem of an extrapolation of the finite-size exact- 
diagonalization results to the larger system sizes. In ID this problem is especially im- 
portant since the finite-size cross-over regions in the vicinities of the critical points have 
only logarithmical macroscopic smallness. (For a typical example see the spectrum in Fig. 
3 demonstrating quite a smooth behavior at V = 2.) Obviously, the extrapolation can be 
performed with the RG equations. As the parameter K naturally follows from the spectrum, 
the only problem is to obtain y. This can be done by fitting the behavior of K as a function 



of N a with RG equations. The details and results of such an approach see in [24 



V. POSSIBLE REALIZATIONS OF ID SUPERCURRENT STATES 

In this section we discuss in some detail different systems and experimental setups which 
allow to study SC states in a ring geometry. The most obvious choice is a real ID electron 
wire of mesoscopic size. Recent experimental techniques p6| produce rings with L ~ 10 4 A° 



and only few transverse levels (in ideal case one) below the chemical potential \i. The 



experiments can be done on ensembles as well as on single rings [27]. In a static experiment 



one may study SC by measuring the so-called persistent current, which is defined as 

j = -^dF/d<p (y? = 0/0o), (57) 

where F is the free energy. At zero temperature j = 2vi/L(I — (p), where I minimizes the 
potential energy Ei(<p). In ideal system j((p) is a saw-like function with abrupt change of 
the current from —vj/L to Vj/L at tp = m + 1/2, when the two energy levels 1 = and 
1=1 are degenerate. The singularity is changed to crossover when we incorporate the finite 
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energy splitting Aab into the problem. A trivial solution of the two-state Hamiltonian gives 
now 

m = g% ( , 2nVl/L - l) ; 6<p = <p- 1/2. (58) 

We see that the crossover region is very narrow | Sip |~ LAab/2vj in the superfluid phase. 
Clearly, such a ring may work as a very sensitive magnetometer capable of measuring the 
flux with the accuracy ~ 4>o(vi/ LAab) <C 4>o- (An interesting experimental setup was 



proposed recently in p8fl , where the ring conductance in the flux is probed by two weak 
tunnel junctions at a distance L/2. In this case the Josephson critical current at zero 
voltage changes dramatically at <p = 1/2.) 

The static picture just described is valid provided the time of experiment is much 
longer than the relaxation time r (see Eq. fl53|) ) , or if the flux was switched on adiabati- 
cally d In ip / dt <C A^. For large enough systems the both requirements might be violated. 
By inserting nonadiabatically the flux <p « m (same arguments go through for ip = 1/2- 
integer) we prepare the SC state \ I = m). Then, according to the standard two-state 
dynamics, the current starts oscillating with the frequency [(8iinvj/ L) 2 + A^] 1 / 2 . In a more 
sophisticated experiment one may study the echo signal from the ring. Actually, the echo 
or the burned hole experiment on the ensemble of rings with different sizes seems to be 
the most appropriate one, because the wide spread in A/ and if for individual rings makes 
it difficult to see the resonant absorption of the electro- magnetic field. We are not going 
into details here because this case is identical (with clear modifications for the distribution 
function) to that of the two-level systems in amorphous materials (see e.g. the review article 



The other system described as a ID Luttinger Liquid at very low temperatures is in fact 
a 2D electron gas in the FQHE state. In the FQHE state all bulk excitations are frozen 
at temperatures much less than the energy gap in the collective spectrum of interacting 
electrons, and the only low-energy degrees of freedom left in the incompressible electron 
liquid are the edge excitations which form ID chiral Luttinger Liquid f30H33f. Here we will 
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discuss the sample with the Carbino disk geometry where the right- and left-moving fields are 
described by the inner and outer edge currents correspondingly, and the relaxation process 
is via quasiparticle tunneling between the edges through the bulk at the point of constriction 
(or many constrictions if we are interested in the disordered case). The crucial parameter 
K~ l in this case equals to the filling factor of the FQHE state v [[5tJ, and for v = l/(2k + 1) 
(when there is only one channel on each edge) we are in the insulator phase, that is in the 
macroscopic sample L — > oo the constriction is equivalent to cutting the ring. One may 
argue that our results do not apply to this case because we concentrated on the superfluid 
phase only. However in all our calculations done in sections |J, |TXT] the only assumption used 
was that the dimensionless vertex w is small. In the superfluid phase this was true even for a 
large microscopic " bare" vertex. Suppose now that we have a finite-size ring with extremely 
small "bare" value of w, so that even for large K the scaling toward larger values of w still 
leaves w(L) small, and we actually never enter the true insulating phase. In this case all 
our results remain intact and directly apply to the experiment. Moreover, large values of 
K might be an advantage because we need a ring with the distance between the edges d 
much large than the magnetic length l H , in order to be able to speak about well-defined 
edge excitations separated by inert bulk. The tunneling rate w goes to zero exponentially 
with the distance d, and for not too narrow disks will be unobservable if not for the increase 
of w under the scaling when K is large. 

The other advantage of the Hall system is that it is not sensitive to small deviations of 
the magnetic field from special values, and the width of the disk d does not constitute any 
problem when an external magnetic field is varied to change the flux through the inner edge. 
Obviously, the above discussion of possible experimental setups for studying the dynamics 
of supercurrent states is in order here. One can also make use of the spatial separation 
between the left- and right-moving currents (which is not the case in really ID wires). In 
the Hall system current relaxation is connected with the charge transfer between the edges, 
which allows to follow, e.g., current oscillations, by monitoring the charging effects in the 
vicinity of the sample. 
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In all examples discussed so far the parameter K was larger than unity (in ID conducting 
wires due to Coulomb repulsion between the electrons we expect K > 1). Consider now a 
3D superconducting ring having the two spatial dimensions d\ and o?2 (defining the sample 
cross-section) much smaller than the ring circumference L. The question we address here 
is the size dependence of the supercurrent relaxation/oscillation in such a ring. The niave 
answer is: since the relaxation is via vortex line tunneling across the constriction (or vortex 
ring nucleation inside the ring with subsequent expansion to the sample boundaries), then 
the process will depend on d\ and d,2, but not on L because it happens locally right at 
the point of deliberately introduced constriction and none of the vortex sizes exceeds the 
cross-section diameter. As we demonstrate below, this answer is absolutely wrong, and the 
relaxation process depends explicitly on the ring length L. 

At temperature T <T C (we assume that the superconducting energy gap is nonzero and 
of order T c all over the Fermi surface) we describe the superconductor by effective action 
for the phase of the order parameter which has the same form as Eq.(^)), except that now 
the integration over dx ought to be understood as dx\ dx<i dx%. Let us go to even lower 
temperatures T <C 1/md 2 (we assume that the cross-section diameter d is larger than the 
coherence length £) , when the transverse excitation modes are frozen, and the effective action 
can be further integrated over dx\ and dx2- What is left is the familiar ID effective action 
with yj 'A s / ' k ~ Vp and A s , k ~ (d/a) 2 . Clearly, in this case we deal with extremely small 
values of K 

K rsj {a/df < 1, (59) 

and the effective ID system is indeed in the superfluid phase. Thus we find that the rate of 
current oscillations will be proportional to 

In A AB /E AB (d/a) 2 In L/d (60) 

according to our general arguments which lead to fl25|) (we used the fact that 1/K ^> 1). 

Till now we completely ignored the effects of Coulomb interaction between the particles, 
so the result (|60| ) rather applies to the neutral superfluid like He-II in a narrow long tube. 
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One might also think about the experiment where the superconducting ring is placed on 
top of the insulator-metal plate, which will screen long-range Coulomb forces. However if 
the interaction between the charges remains unscreened, the result (|B0D has to be essentially 
modified. The spectrum of collective modes in a 3D superconductor in the limit q — > is 
given by ]33J (see also H) 

w 2 ~ v 2 F q 2 U(0)V c (q) ; (61) 

and for Vc(q) ~ e 2 /q 2 has a plasmon gap w 2 = ne 2 /m ~ VpK 2 , where 11(0) is the static 
polarizability, and k" 1 is the normal- metal screening radius. 

Consider now the limit L^ 1 < q <C d' 1 , which corresponds to our quasi- ID geometry 
as discussed above. In this momentum range we almost recover back the sound-like energy 
spectrum for phase fluctuations 

w 2 q ~ v 2 F q 2 (K p d) 2 \nqd . (62) 

Weak logarithmic dependence of the sound velocity and parameter K which is now given as 

K- 1 ~ — —77 (63) 

S« 2 ln 1/2 (qd) K ' 

can not change the theory in any essential way, except that now the integration of the RG 
equation jH^ K~ x dq/q gives In 1 / 2 L/d instead of conventional logarithmic dependence. Thus 
we find for the charged superconductor another law 

]nA AB /E AB — 2 In 1 / 2 L/d. (64) 

Of course, to observe these dependences experimentally the ratio d/a must be not too 
large. Bearing in mind that we assumed d > £, we suggest superconductors with short 
coherence length (e.g., HTSC), or else 4 He-II, although we are not aware of any experimental 
technique for preparing such narrow (d ~ 20 A) long tubes for the helium experiment. 

In this paper we mostly discussed the properties of SC in a sample with the ring topology 
threaded by the flux. We realize, however, that this situation can be modeled in a ID wire 
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connected to the constant-current source through the Josephson contacts. In a way, here 
one deals with a reverse problem of finding the phase difference between the contacts as a 
function of current. This experimental setup might be of particular interest for studing the 
inelastic relaxation of SC due to their coupling to the substrate at large values of the current 
j > jci where j c is the maximum value of the persistent current in Eq.(^). At mesoscopic 
values of the current j ~ the phase slippage process in a ID system is possible only if 
assisted by the energy dissipation to the substrate because we are dealing with the lowest 
energy levels, and ID phonons form a discrete spectrum which is not in resonance with SC 
in the general case. 



VI. ACKNOWLEDGMENTS 

We would like to thank Yu.Kagan and I. Affleck for the discussions and valuable com- 
ments. 

This work was supported by International Science Foundation (Grant No. MAA300) 
(NVP and BVS) and by Russian Foundation for Basic Research (Grants No. 94-02-05755a 
(VAK and AIP) and No. 95-02-06191a (NVP and BVS)), and partially by Grants No. 
INTAS-93-2834 [of the European Community] (BVS) and No. NWO-07-30-002 [of the Dutch 
Organization for Scientific Research] (NVP and BVS). 



33 



FIGURES 



Fig.l. Spectrum evolution of a disordered system with increasing disorder. N a = 11, 
number of particles N b — 7, U — 2.5. Here (2, 0) and (0, 2) are the states with two identical 
minimal-momentum phonons; (1,1) is the state with two minimal-momentum phonons, 
moving in the opposite directions. 

Fig. 2. Spectrum evolution of a disordered system with increasing interaction. N a = 11, 
N b = 7, W = 0.5. The notation is the same as in Fig. I. 

Fig. 3. Spectrum of the commensurate hard-core model as a function of the nearest-neighbor 
interaction; N a = 22, N b = N a /2. 

Fig. 4. Scaling of relative AB-splitting with Na in the case of single impurity (ei = 0.5); 
N a = 8, 10, 12, 14, 16, 18; N b = N a /2. 

Fig. 5. Scaling of relative AB-splitting with Na in the disordered case (W = 0.5); N a = 
8, 10, 12, 14, 16, 18; N b = N a /2. Each point represents an average over 100 realizations. 

Fig. 6. Scaling of the relative splitting of the first SC with Na in the commensurate case; 
Na = 10, 12, 14, 16, 18, 20, 22; N b = N a /2. 

Fig. 7. Scaling of the relative splitting of the first SC with Na in the commensurate case at 
the critical point. iV a = 10, 12, 14, 16, 18, 20, 22, N b = N a /2. 

Fig. 8. The ratio R between the first SC energy splitting A 1 and the square of the lowest 
splitting A AB ; N a = 11, N b = 7, and W = 0.25. 
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